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We report a way of obtaining a spin configuration snapshot, which is one of the representative 
spin configurations in canonical ensemble, in a finite area of infinite size two-dimensional (2D) 
classical lattice models. The corner transfer matrix renormalization group (CTMRG), a variant 
of the density matrix renormalization group (DMRG), is used for the numerical calculation. The 
matrix product structure of the variational state in CTMRG makes it possible to stochastically 
fix spins each by each according to the conditional probability with respect to its environment. 
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§1. Introduction 

For a classical lattice spin system that interacts with a 
reservoir, the spin configuration observed at an instant, 
which we call snapshot in the following, is one of the 
representative configurations in the canonical ensemble. 
Such a snapshot is experimentally observed if the time 
scale of the observation is much shorter than that of the 
evolution of the system. A freezed spin configuration 
after sudden cooling can also be regarded as a kind of 
snapshot. In general, snapshots show rough outlook of 
the system. For example, the typical size of a spin in- 
verted island in the ordered state is of the order of the 
correlation length, and symmetries of ordered states can 
be intuitively identified. Figure 1 shows a snapshot in- 
side the area of 100 by 100 of the two-dimensional (2D) 
ferromagnetic Ising model, where the system size is much 
larger than the shown region. Normally such a snapshot 
is drawn by Monte Carlo simulation, ^ while the shown 
one is created by the corner transfer matrix renormaliza- 
tion group 7 - 1 (CTMRG), a variant of the density matrix 
renormalization group 2 ) (DMRG) applied to 2D classical 
lattice models. So far both DMRG and CTMRG have 
been used for calculations of thermal average of spin cor- 
relation functions, but not for the snapshots. 

In this article we report a way of obtaining the con- 
ditional probability for a row of spins of length N sur- 
rounded by the rest of the system, using the matrix prod- 
uct structure 3 ) (MPS) of the variational state employed 
in CTMRG and DMRG. These spins can be fixed one 
by one according to the conditional probability, and af- 
ter fixing all the spins in the row it is possible to obtain 
the conditional probability for the next row in the same 
manner. Applying such a fixing process for M numbers 
of rows, one obtains a snapshot for the area of N by M 
in infinite — or sufficiently large — 2D lattice systems. 
Numerical cost for this snapshot observation is of the 
same order of conventional zipping process in the finite 
system DMRG algorithm. 2 ' 4 ) 



Fig. 1. Snapshots of the ferromagnetic Ising model at T = 1.5 
(left) and 2.27 (right), where we have chosen the nearest neighbor 
coupling constant as the unit of energy. 



in the 2D lattice model in terms of the corner transfer 
matrix formalism. In section 3 we generalize the spin 
fixing procedure to a row of spins of the length N, by 
taking partial sum for the inner product between MPSs. 
Introducing position dependence to the local factors that 
construct MPS, we successively fix the M-rows of spins 
as we explain in Sec. 4. In the last section we conclude the 
result and discuss the relation with quantum observation 
in one dimension. 

§2. One Spin Fixing 

The corner transfer matrix (CTM), which was in- 
vented by Baxter, 5 ) is not only useful for rigorous anal- 
yses of 2D lattice systems, but also efficient for numer- 
ical calculations of thermodynamic functions, especially 
away from the critical point. Let us briefly look at the 
construction of CTM and the block spin transformation 
applied for it. 

We consider the square lattice Ising model as an exam- 
ple of the 2D classical lattice spin models. The partition 
function of the system is expressed as 



Z = ^exp(ffM), 



(2.1) 



In the next section we briefly explain how to fix a spin 



2 



Kouji Ueda, Ryota Otani, Yukinobu Nishio, Andrej Gendiar, and Tomotoshi Nishino 



C ^ B 



! d $ a : 



a 



(a) 



1 c 


B i 


H i 1 v 


, D 


a ; 





(b) 



Fig. 2. Division of the 2D lattice into 4 corners, (a) Variables of 
the corner transfer matrix A. (b) Rcnormalized corner transfer 
matrices. 



where {a} represents a configuration of Ising spins a = 
±1 on the square lattice, and H{a} is the Ising Hamil- 
tonian, that is a sum of Ising interaction over all the 
bonds. 

In the CTM formalism, the 2D lattice is divided into 
4 parts, so called the quadrants or the corners. 5 ' Let us 
label each quadrant as A, B, C, and D. The Hamilto- 
nian H{a} is then expressed as sum of the corresponding 
parts 

H A {a A } + H B {a B } + H c {a c } + H D {a D } , (2.2) 

where {cr A }, {cr B }, {i c }, and {a D } denote spin configu- 
rations in each quadrant of the system. Note that neigh- 
boring quadrants share the same spins at the boundary 
between them. In this way the Boltzmann weight of the 
whole system is expressed as a product of 4 factors 

exp(#{cr}) = exp(H A {a A }) exp(H B {a B }) 

exp(H c {a c }) exp(H D {a D }) . (2.3) 

The corner transfer matrix is the partial sum of each 
factor with respect to the spins inside the corner 



A(aa(3) = ^ exp(H A {a A }) 



(2.4) 



except for those spins at the boundary with other cor- 
ners, where the notation £' denotes this restricted sum- 
mation, and where a and (3 denotes groups of spins at 
the boundary. (See Fig. 2 (a).) The CTM is a block diag- 
onal matrix with respect to a, while in this article we do 
not use this matrix property explicitly. Other CTMs B, 
C, and D are defined in the same manner. If the system 
is invariant under the rotation of 90 degree these CTMs 
are equivalent. For simplicity, we assume this symmetry 
in the following. 

The matrix dimension of the CTM increases exponen- 
tially with the system size. In order to avoid this blow up 
and treat the CTM accurately in numerical calculation, 
block spin transformation is introduced in CTMRG from 
those spins a and (3 at the boundary to m-state auxiliary 
variables. Such a transformation maps the CTM A into 
the compressed (or renormalized) one 

MS™) = E W U M A ( aa ® > ( 2 - 5 ) 

a/3 



where £ and v are the m-state auxiliary variables. The 
transformation matrix (a) is obtained by diagonal- 
izing the density matrix p = ABCD. 2 '^ Other CTMs 
can also be mapped to B(ucrrj) 7 C(r]a^i), and D(fia^) as 
shown in Fig. 2 (b). The approximate partition function 



Z' = J2 A (&v) B{vari) C(r]afi) D(na^) 



(2.6) 



is close enough to the original one Z in Eq. (2.1) if m 
is sufficiently large; normally m is of the order of 10 to 
1000. 

When the system size is far larger than the correlation 
length of the system, the renormalized CTMs becomes 
independent of the system size except for a constant mul- 
tiple. Throughout this article we assume such a condi- 
tion, and regard the renormalized CTM as a quadrant of 
infinite size systems. 
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Fig. 3. The left and the right vectors created from CTMs. 

Now we calculate the probability of observing a = 1 
(up) and a = — 1 (down) at the center of the system. We 
first prepare two partial sums 

= £>(£»/) 5(1/^), ( 2 - 7 ) 

that we regard 2m 2 -dimensional vectors in the following. 
Figure 3 shows the graphical representation of the above 
equation. The order of variables in L(rja^) is opposite 
to that of R(£<jrj), since we keep the order of auxiliary 
variables in the right hand sides of the above equations. 
The approximate partition function Z' is represented as 

Z' = {L\R) = J2L(v^)m^v), (2-8) 

£<77J 

where we have introduced bracket notations for the book 
keeping. The probability of a taking the specific value a 
can be written as 

(L\S(a,a)\R) _ £ £ „ L(rp£) R^ar]) 



p{a) 



(L\R) 



(2.9) 



where the probability satisfies the normalization p(l) + 
p(— 1) = 1. Therefore, creating a random number x in 
the range [0,1) one can fix a to a = 1 if x < p(l), 
otherwise to a = — 1. 
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§3. Snapshot in a row 

Let us generalize the spin fixing procedure to a N num- 
bers of spins in a row from a x to a N , where TV is much 
smaller than the system size. For this purpose we intro- 
duce the half- row transfer matrices (HRTMs), that are 
upper and lower halves of the transfer matrix T to the 
horizontal direction. Figure 4 shows the system that we 
consider in this section. Between the vectors L and R, 
there arc N — 1 numbers of transfer matrix 



T 



Sa Pa 



(3.1) 



from i = 1 to TV - 1, where S t and P i are HRTMs 

Si = S(£i <T t <7 i+1 £ i+1 ) 

Pi = p (Vi+i^i+i^iVi)- (3-2) 

We have aligned the variables of HRTMs in the clock- 
wise order as CTMs. At the moment all the HRTMs 
are equivalent, though we put indices for identification 
of their positions. 




R 



& & & 



Fig. 4. Half-row transfer matrices between vectors L and _R. 



We fix these N spins from the left to the right, 
successively calculating the conditional probability 
Ps 1 ...s. {&i ) f° r the spin at i-site after fixing those spins 
in the left a 1 , . . . ,a i _ 1 . The way of calculation is es- 
sentially the same as operator multiplication to a given 
matrix product state 3 ~ 5) (MPS) in finite system DMRG 
algorithm, and it is numerically important to prepare 
partial sum of local factors in advance. We write R as 
R N since it contains a N as its variable. Let us multiply 
the transfer matrices T i one by one to the vector R N . 
We obtain 



Rn-i — ^JV-l 

Rn-2 = Pn-i i etc., 



(3.3) 



down to R 1 . In the numerical calculation we do not 
possess the transfer matrix T i explicitly, but we apply 
S i and P i part by part as 



R(ViViti)= s (£i a i a i+i£i+i) 



(3.4) 



Y R{^i+i^i+iVi+i)P(Vi+i^i+i^iVi) , 

where the second sum should be taken first. (See Fig. 5.) 
After obtaining R 1 we can calculate the probability for 



L Oj >k— 6 a 2 ct n _, 6—4 ct ;V r 



Fig. 5. Multiplication of Transfer Matrix. 



a x as before 



(3.5) 



(iil^i) 

where we have written L as L x since it contains a x as 
its variable. According to this probability we stochasti- 
cally fix the first spin in the row. We then consider the 
conditional probability p s (<r 2 ) for the second spin after 
fixing the first one. This time, we have to prepare L x 
multiplied by T 

(3-6) 

as shown in Fig. 5. It should be noted that we do not 
take configuration sum for a x since it is already fixed; L 2 
contains a fixed spin a x in it. Using R 2 calculated before 
and L 2 in Eq.(3.6), we obtain the conditional probability 
for the second spin 



(L 2 \5(a 2l a 2 ) \R 2 ) 



(3.7) 



(L 2 \R 2 ) 

According to this probability we can fix the second 
spin. After that we calculate L 3 in the same manner 
as Eq.(3.6). Repeating these spin fixing procedure to 
a N , we finally obtain a spin snapshot {a} = a x , cr 2 ,..., 
a N for a group of N spins in a row. 

§4. Snapshot in a Rectangular Area 
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Fig. 6. The double-row spin system. The cross marks represent 
the already fixed spins in the first row, and the circles are those 
spins that we will fix each by each. 



Let us consider the way of fixing M numbers of spin 
rows successively. The first step is to obtain the proba- 
bility for the second spin row under the condition that 
the first row a x , a 2 , ■ ■ ., <J N is already fixed. For the dis- 
tinction let us write the spins in the second row as {t} = 



4 



Kouji Ueda, Ryota Otani, Yukinobu NlSHIO, Andrej Gendiar, and Tomotoshi NlSHlNO 



r l5 r 2 , . . ., Tjy. Figure 6 shows the system that we con- 
sider for a while, where there is a transfer matrix to the 
vertical direction between {a} shown by cross marks and 
{r} by circles. This transfer matrix consists of the right 
HRTM 0(ya 



N N 



the left one Q(fj,'T 1 o 1 n), and N — 1 
numbers of local Boltzmann weights W{T i T i+l d i o i+ i ) 
from i = 1 to N — 1 in between. 8 ) 

We reduce the above double-row system to the single- 
row one treated in the previous section by way of exten- 
sion of the CTMs and HRTMs to the vertical direction, 
similar to the system size extension in CTMRG. 7 ) The 
HRTM S is extended by putting W on top of it 

= W(r l T l+1 a l a l+1 )S(£ l o l o l+1 £ l+1 ) 

(4-1) 

as shown in Fig. 7. Here we do not regard a i and <J i+1 
as the variable of the extended HRTM S' since these are 
fixed constants. Thus the number of elements in S[ is the 
same as that of S i . It should be noted that after such 
an extension S[ becomes position dependent. 
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Fig. 7. Extension of the HRTM S to the vertical direction. 



The area extension of CTMs A and D are done by 
putting HRTMs O and Q 



a '{£,nTn v ') = A ^n°n v ) 0{vo n t n v') 



£>'(Ai£i) = E Difio^) , (4.2) 

respectively, as shown in Fig. 8. Again, the number of 
elements in the extended CTMs are the same as the orig- 
inal ones. In this way we have reduced the double row 
system in Fig. 6 as single row one constructed from A', 
B, C, D', P i , and S 1 -. The spin fixing procedure from 
the left to right explained in the previous section is now 
applicable to the second row of spins {r}. It is straight 
forward to repeat the extension of HRTMs and CTMs to 
the vertical direction for M — 1 times, we finally obtain 
the snapshot of the size N x M in the infinite system. 
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Fig. 8. Area extension of CTMs D and A. 



snapshot observation can be applied to various 2D lat- 
tice models with short range interaction, such as the IRF 
and the vertex models. 9 ) It is also possible to treat half- 
infinite or finite size systems and observe snapshot near 
the system boundary if we admit position dependence to 
HRTMs. (The finite system DMRG may be more appro- 
priate than CTMRG in such a position dependent case.) 

The calculations of conditional probability requires 
renormalized CTM and HRTM, that are converged to 
the infinite system size limit in the numerical algorithm 
of CTMRG. The convergence is quite slow near the criti- 
cal temperature, and it is necessary to accelerate it. The 
same problem exists in the infinite system DMRG algo- 
rithm, and the problem has not been solved yet. 

It is possible to apply the conventional real space 
renormalization group transformation, such as the ma- 
jority rule for blocks of spins, to the obtained snapshots 
and calculate the critical indices. This suggests that both 
CTM and HRTM implicitly have information about crit- 
icality such as critical indices and scaling functions. To 
pull out these directly from CTM and HRTM, without 
looking at the snapshot is one of the future problem. 

We finally comment about a relation between 2D clas- 
sical systems and ID quantum systems. The snapshot in 
the former corresponds to many time observation result 
in the latter. A snapshot can be regarded as a repre- 
sentative path in the path integral formulation of quan- 
tum systems. Successive observations in time-dependent 
DMRG 10,11 ) for numbers of time slices is the same type 
of computation. 
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§5. Conclusion and Discussion 

We have explained the way of obtaining spin snap- 
shot in the area of N x M for the 2D Ising model. The 



